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Abstract Stochastic mean-field simulations for multifragmenting sources at the same excitation energy per 
nucleon have been performed. The freeze-out volume, a concept which needs to be precisely defined in this 
dynamical approach, was shown to increase as a function of three parameters: freeze-out instant, fragment 
multiplicity and system size. 

PACS. 2 1.30.Fe;25.70.-z;25.70.Lm;25.70.Pq;24.60Ky; NUCLEAR REACTIONS 119 Sn( 129 Xe, X), E = 32 
AMeV; 238 U( 155 Gd, X), E = 36 AMeV; central collisions; multifragmenting sources; stochastic mean-field 
simulations; dynamical evolution of the freeze-out volume. 



1 Introduction 

Several hundreds of nucleons may be brought into inter- 
action in central heavy ion collisions around the Fermi 
energy. Such reactions are good candidates to provoke a 
liquid-gas type phase transition - conceivable given the 
nature of the nuclear force - and to break the system into 
massive fragments QH]. The volume of such a source of 
ejectiles at the instant when all of them become free of 
the attractive force and feel only the Coulomb repulsion - 
the freeze — out volume - brings information on the coex- 
istence of phases. It is a key quantity |Sj to be connected 
to the physical observables, asymptotically measured. If 
in the statistical models (31011121 it is a basic apriori hy- 
pothesis, in the following dynamical treatment this volume 
is provided - at a given available energy - as a family of 
results illustrating the temporal and spatial evolution of 
the source in multifragmentation. 

Nuclear multifragmentation may occur when the source 
has expanded through the spinodal region of negative com- 
pressibility pQ of the liquid-gas coexistence domain, a sce- 
nario valid for other many-body systems too 8,9. Related 
Stochastic Mean Field (SMF) approaches [1(111111112) con- 
sider, under different approximations, the amplification of 
density fluctuations, due to N-body correlations, by the 
unstable mean field leading to spinodal decomposition. 
The Brownian One-Body (BOB) dynamics version, simu- 
lating the fluctuations by means of a brownian force in the 
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mean field [13114115) . coupled to Boltzmann-Nordheim- 
Vlasov (BNV) one-body density calculations ^3] , was suc- 
cessfully confronted with multifragmentation data mea- 
sured with the 47T multidetector INDRA ^H] . Two systems 
at close available energy per nucleon ~ 7 MeV, were exper- 
imentally studied: 32 AMeV 129 Xe + na 'Sn and 36 AMeV 

tween simulated events for central collisions - filtered ac- 
cording to the experimental features of INDRA - and the 
experimental data was recently extended from fragment 
multiplicity M, charge Z, largest charge Z max and bound 
charge Zb om d = SZ distributions ^BJ to charge [TT?ll2T] 
and velocity correlation functions and energy spectra |22j . 
The theoretical approaches developed in relation with the 
advanced experimental methods of nuclear physics may 
be esteemed in the new and more general physics of the 
phase transition in Unite systems |3ll51l23) . 



As a reasonably successful dynamical description of 
the multifragmentation process at intermediate bombard- 
ing energies, the above mentioned three dimensional (3D) 
SMF simulations provide a well adapted framework to 
address the question of the freeze-out volume. We focus 
in the present paper on disentangling the time, fragment 
multiplicity and source size of the freeze-out volume de- 
pendence. A pragmatical definition of the freeze-out in- 
stant is proposed. 
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2 Expanding sources provided by BOB 
simulations 



SMF simulations of nucleus-nucleus collisions based on the 
Boltzmann-Langevin equation were proposed to treat un- 
stable systems [TTIll 21124] but their application to 3D nu- 
clear collisions is prohibited by computational limitations. 
Instead, the spinodal decomposition of two diluted nuclear 
sources, mentioned above, at the same temperature (~ 4 
MeV), was mimicked through the BOB dynamics 1 1311 151 
IT8| . applicable to locally equilibrated systems. The brow- 
nian force employed in the kinetic equations is grafted on 
to the one-body density evolution, calculated in a BNV 
approach |14| . at the moment of maximum compression 
to ~ 40 fm/c after the beginning of the collision, before 
the entrance of the system into the spinodal region. Its 
strength can be tuned to describe the growth of the most 
important unstable modes, ascribed to the density fluc- 
tuations, which need a short, but finite time to develop. 
The dispersion relation |25| puts them in evidence. It in- 
cludes quantal effects and connects the characteristic time 
to its associated multipolarity. In turn, the multipolarity 
of the unstable collective modes, increasing with the size 
of the source, may be related to the fragment multiplicity 
[251126) . The delimitation between fragments - the "liquid 
drops" - and light clusters - "the gas" - in which they are 
embedded is given by a cut-off value p m in > 0.01/m -3 of 
the nuclear density p |18| . 

The ingredients of the simulations, corresponding to 
zero impact parameter, as well as the selection criteria 
for the complete events from experimental central colli- 
sions, can be found in |18II27| . The fragments are defined 
as having the atomic number Z > 5. In the reported cal- 
culation, as in the experimental selection, only events hav- 
ing the final fragment multiplicity M > 3 were considered 
|18| . The calculated total charge of the multifragment- 
ing sources in the spinodal zone, Z tot = 100 for 129 Xe + 
119 Sn and Z tot = 142 for 155 Gd + 238 U are close to the 
experimentally reconstructed ones ^Hj; INDRA identifies 
the mass of the light charged products [Z < 5) but not 
that of the fragments. The Skyrme force used in our sim- 
ulations is not isospin dependent. Consequently, the mass 
numbers of the sources: A to t — 238 (for the light one) and 
Atot = 360 (for the heavy one) correspond to the con- 
servation of the entrance channel N/Z ratio. The charge 
distributions (normalized to the fragment multiplicity of 
the event) measured j2H] or simulated in BOB calculations 
|18| . are identical in the two cases - consistent with a bulk 
effect in the involved multifragmentation process. 

Starting from the two initial partners of the reaction, 
the BNV calculation leads to a spherical distribution of 
matter. It undergoes a self-similar expansion, which is 
not dramatically altered by the BOB simulations - Fig. 
^ The simulated sources continue their rather isotropi- 
cal 3D expansion in time, no particularly elongated or flat 
shapes being produced. The concept of a radially symmet- 
ric freeze-out volume keeps its full sense, provided that the 
related instant may be singularized. 
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Figure 1. One event density evolution for each of the two 
multifragmenting sources: A to t = 238 - upper panel - and A to t 
= 360 - lower panel. The collision direction in the entrance 
channel is along Z axis and the unit is fm on both axes of 
these X - Z views. The centre of mass coordinates are X = Z 
= 20 fm. 

2.1 The freeze-out instant definition 

The fragments are not all formed at the same moment, 
their mean multiplicity increasing in time, up to ~ 250 
fm/c ^H! when it saturates at a value of about 5 for the 
lighter source (Z to t — 100, A to t = 238) associated to the 
Xe + Sn reaction, and of about 8 for the heavier one {Ztot 
= 142, Atot = 360) associated to the Gd + U reaction. 
Indeed, even for the same final multiplicity of a source, 
there are events where the density fluctuations grow up 
faster and others where the process is slower. A question 
of definition appears related to the freeze-out instant. Our 
BOB calculations are recorded in steps of time At = 20 
fm/c starting from to. Each event is traced back in steps of 
20 fm/c from its asymptotical configuration (t w 250 fm/c) 
of final multiplicity AI to the moment when the fragment 
multiplicity decreases one unit. We define the freeze-out 
instant of an event as the moment when its final fragment 
multiplicity M was established. It means that if at t^i — 
to + (i — l)At the event has the multiplicity M — 1 and 
at ti = to + iAt it reaches the final multiplicity M, ti will 
be considered as freeze-out instant. All events which got 
their final multiplicity M at the moment U are treated 
together. Fig. [5] shows, as an example, the distribution 
of the moments t = ti at which the final multiplicity M 
= 5 was reached for the lighter source A to t = 238: t e 
[120,260] fm/c. 

2.2 Freeze-out configurations 

For the sake of simplicity, the fragments are considered 
as spheres of normal density, but the results are quite 
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Figure 2. The distribution of the moments at which the final 
multiplicity M = 5 was reached for the source having Atot = 
238. 



independent of this particular hypothesis. The relative 
distance djk between the surfaces of two fragments j, k 
(j,k = 1t¥ and j ^ k) in one event, that will be called 
the intra-fragment distance in the following, is much vary- 
ing for a given final multiplicity M: from a minimum of 
the order of 1 fm between the two most recently split- 
ted fragments to a maximum value which increases with 
t. This behaviour is shown - for the source (100,238) and 
the multiplicity M = 5 - in the left-upper and left-middle 
panels of Fig. |3| at two different freeze-out instants: 180 
and 240 fm/c. The distribution evolves towards larger dis- 
tance values in time, from an asymmetric to a more sym- 
metric shape. For the same two instants, the right-upper 
and right-middle panels of Fig. present the distributions 
concerning the source (142,360) and M = 8. The evolu- 
tion is similar to the lighter source case. Together with the 
lower panels, related to the freeze-out instant t = 180 fm/c 
and the multiplicity M = 7 for both sources, the upper 
graphs let also one see that, at the same moment (t = 180 
fm/c), the mean value of the relative distance depends on 
M for a given source, while for different sources and the 
same multiplicity, e.g. M = 7, on the size of the source. 

The distribution of the intra-fragment distances at a 
certain moment is quite illustrative for the momentary 
spatial configuration of the source: the shortest distances 
concern the first order neighbours, the longest ones the 
largest size of the sources and, above all, the intermediate 
values are qualitatively informing about the dilution of the 
source. The first and second order moments of the distri- 
bution and their dynamical evolution are appropriate to 
synthetize this last aspect. As expected, these quantities, 
analysed multiplicity by multiplicity, are increasing with 
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Figure 3. The distributions of relative distance djk between 
the surfaces of two fragments j, k, simulated event by event 
and related to two sources: Atot = 238 - left column - and 
Atot = 360 - right column, for various freeze-out instants and 
multiplicities. 



increasing t = tj, testifying on the enlargement of the mat- 
ter distribution in the nuclear source in time. These values 
are higher for the heavier source than for the lighter one, 
and their increase in time seems to be slightly more pro- 
nounced. The fact has to be put in connection with the 
total Coulomb repulsion and the radial flow, higher for the 
heavier system than for the lighter one |18U22| . An inter- 
esting result is that, at constant values of t = ti, the mean 
and the variance of intra-fragment distance distributions 
increase with the multiplicity - a little bit more accentu- 
ated for the lighter system. The slope of this dependence 
is higher at later moments t. 

For each group of events reaching a final multiplicity 
M at a certain moment t = t», one may look for the corre- 
sponding local fragment concentration: <1N /&V as a func- 
tion of the vector radius absolute value r of the fragment 
position in the source reference framework. Examples of 
such distributions (at t = 180 and 240 fm/c) are given 
in the left-upper and left-middle panels of Fig. 01 for the 
lighter source and M — 5. When the same source and final 
multiplicity are involved, the mean radius f and full width 
at half maximum (FWHM) increase with the time t = U 
when the last two fragments are separated. In the right- 
upper and right-middle panels of Fig. 0] are represented - 
for the heavier source - the local fragment concentration 
distributions at these two moments and the final multi- 
plicity M = 8. As previously, the distribution evolves in 
time towards larger radii; these radii are longer than in 
case of the lighter source. For the same final multiplicity: 
M = 7 and the same moment: t = 180 fm/c - the distri- 
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Figure 4. Radial distributions of the fragment concentration, 
simulated event by event and related to two sources: A to t = 
238 - left column - and Atot = 360 - right column, for various 
freeze-out instants and multiplicities. 



butions represented in the lower panels for both systems 
- the mean radius f and FWHM are longer in case of 
the heavier source. The upper and the lower pannels to- 
gether show that, for a given source at the same freeze-out 
instant, f varies in the same sense as M. 

Except the lowest multiplicity case for the lighter source, 
the radial distributions of the fragment concentration at 
the freeze-out instant are practically empty towards short 
radii, as shown in Fig. fj] They may be interpreted hence 
as reminiscences of bubble-like configurations. 
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Figure 5. Freeze-out volumes as a function of the freeze-out 
instant, for two sources: A to t = 238 - left column - and A to t = 
360 - right column. The full curves concern full spheres, while 
the dashed ones correspond to hollow spherical envelopes (see 
the text for explanations). 



3 The freeze-out volume 

Once the freeze-out instant defined and the configurations 
when the latest formed fragments get free of nuclear inter- 
action found, one may look for the corresponding volume. 
A sphere of radius f + FWHM/2 englobes most of the 
fragment centres. Its volume - a good estimate of the quan- 
tity of interest - normalized at the volume Vq = (1-2) 3 A to t 
fm 3 of the source at normal density po, is considered as 
a function of the freeze-out instant t — t, at constant fi- 
nal multiplicity M of fragments. It is represented with full 
curves in the left column panels of Fig. [3] for the source 
A to t = 238. The expansion of the source in time is ev- 
idenced: the dynamic process delivering the same final 
number M of fragments implies higher volumes if it takes 
place at later instants. These volumes increase from top to 
bottom with M, as well as the slopes of their variation. By 
adding to the above spheres the complements of fragment 
volumes which are exceeding them, one gains roughly 10% 



on the freeze-out values, but the general behaviour is not 
changed. 

In fact, as shown in Fig. 01 the shapes of the frag- 
ment concentration distributions at freeze-out are gener- 
ally gaussian like. About 75% of their integral is hence 
comprised between f - FWHM/2 and f + FWHM/2. 
The hollow envelopes delimited by the spheres with these 
radii are a kind of lower limits of the corresponding freeze- 
out volumes; always normalized at Vq, they are drawn as 
dashed curves. The right column panels of Fig. show 
similar results for the source A to t — 360. The curves seem 
to rise slightly more rapidly than in the left column, in 
relation with the Coulomb effect and the radial flow. 

The evolution with the multiplicity of the freeze-out 
volume defined above, calculated for full or hollow spheres 
and normalized at the volume of the corresponding source 
at normal density, is shown in Fig. at given freeze-out 
moments. The dilution of the source increases with the 
fragment multiplicity - slightly faster for the lighter sys- 
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fragments, the larger the number of zones between them. 
Consequently, at the same freeze-out instant, the higher 
the multiplicity, the bigger the source volume. From the 
energetic point of view, a higher final fragment multiplic- 
ity M implies a higher fraction of the excitation energy 
consumed as binding (mainly surface) energy. 

The freeze-out volumes provided by the present BOB 
calculation, averaged over time and multiplicity, are, of 
course, well fitting with those extracted, in the same frame- 
work, by using the mean multiplicity information. The 
corresponding densities p are compatible with the gen- 
eral prediction po/10 < p < po/2 of the Statistical Model 
for Multifragmentation (SMM) 0, in particular with the 
value po/3 used to study - in a nonsphcricity hypothesis in 
this latter model - the fragment velocity correlations in the 
32 AMeV 129 Xe + ™ at Sn system [2HJ. The present density 
values are lower than the average densities corresponding 
to the same domain of excitation energy, extracted from 
nuclear caloric curves |.'<()U.'il j . The microcanonical model, 
analysing - all multiplicities together - the same multifrag- 
menting systems 0, leads to volume values close to the 
present results but with a slightly lower dilution for the 
heavier system than for the lighter one. 

The increase of the freeze-out volume with the basic 
experimental observable which is the fragment multiplicity 
in an event is more pronounced than that obtained when 
compacity criteria are used to fill the freeze-out domain 
in a static scenario. For the first time, reliable values as 
a function of time and final multiplicity are calculated. 
They give a dynamical image of the multifragmentation 
process at Fermi incident energies. 

The final multiplicity of fragments is a measurable ob- 
servable, while the time information is not directly ac- 
cessible in the experiments. By weighting the freeze-out 
volume values got at different freeze-out instants with re- 
spect to the corresponding number of events, one may 
obtain a kind of average freeze-out volume at a given final 
multiplicity. This "time" averaged quantity, which looses 
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Figure 6. Freeze-out volumes as a function of multiplicity, at 
different freeze-out instants, for two sources Atot = 238 - left 
column - and A to t = 360 - right column. The full curves con- 
cern full spheres, while the dashed ones correspond to hollow 
spherical envelopes (see the text for explanations). 




Figure 7. As in Fig. |S| but for "time" averaged freeze-out 
volumes; see the text for explenations. 



tern, presenting a higher relative variation of M (left col- 
umn panels) than for the heavier one (right column pan- 
els); the variation is more pronounced at larger times. It 
reflects the mechanism considered here for the multifrag- 
mentation: the density fluctuations. The separation of the 
fragments in such an expanding source is reached on be- 
half of the lower density domains: the larger the number of 



much of its physical signification, indeed, is rapidly rising 
with M, as shown for the two sources in Fig. [7| The vari- 
ety of fragment configurations, connected to the dynam- 
ical fluctuations of the system volume on the fragmenta- 
tion path, has been recently related to monopole oscilla- 
tions |Hj ■ They may push the system towards a metastable 
configuration which eventually recontracts leading to low 
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fragment multiplicities, or develops into a hollow configu- 
ration which fragments at higher multiplicities. 

4 Conclusions 

Stochatisc mean field simulations concerning the multi- 
fragmenting sources formed in 32 AMeV 129 Xe + 119 Sn 
and 36 AMeV 155 Gd + 238 U reactions have been per- 
formed. Their validity has already been checked [18111911211 
12*2*] by reproducing measurable physical observables deter- 
mined for these reactions with the 4ir multidetector IN- 
DRA. We found that the moment of separation of the 
latest two nascent fragments - the definition adopted for 
the freeze-out instant - is mainly distributed in the range 
~ (150 — 250) fm/c; multifragmentation is a dynamical 
process which is fast but needs a finite time. The topol- 
ogy of the associated freeze-out configurations are more 
complex but also more realistic as compared to the sim- 
plifying hypothesis in which the fragments, separated by 
a distance of the order of the nuclear interaction range, 
are forced to fill in a prescribed volume, as generally done 
in statistical codes. 

The corresponding freeze-out volume could thus be dis- 
entangled in connection with the freeze-out instant, final 
fragment multiplicity and source size. For a given source 
and a given multiplicity, the evolution of this quantity il- 
lustrates the continuous expansion of the source in time. 
On the other hand, at the same freeze-out instant, the vol- 
ume of one source increases with the fragment multiplicity, 
a basic measurable observable. As part of the energy as- 
sociated to the internal degrees of freedom of the source is 
going along into fragment separation energy, its volume is 
dramatically increasing. And finally, for the same moment 
and the same multiplicity, the freeze-out volume is bigger 
for the heavier source, involving higher Coulomb repulsion 
and radial flow. 

The dilution of a multifragmenting source, quantified 
as the ratio between its volume at freeze-out and its vol- 
ume at normal density, is therefore increasing with time, 
fragment multiplicity and source size. Further SMF sim- 
ulations, employing an isospin dependent nuclear force, 
would allow to study the isospin fractionation 32,33 , a 
phenomenon expected to be amplified at high source dilu- 
tion. A comparative experimental investigation of isoscal- 
ing characteristics of these two systems prepared at the 
same available energy is desirable too. 
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